NACA TN 2093 


GOVT. DOC 


X , 1 . , , 


y iM'ii/sit/jafz 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2093 


FORMULAS AND CHARTS FOR THE SUPERSONIC LIFT AND DRAG OF 
FLAT SWEPT-BACK WINGS WITH INTERACTING 
LEADING AND TRAILING EDGES 
By Doris Cohen 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 



Washington 
May 1950 




NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2093 

FORMULAS AND CHARTS FOR THE SUPERSONIC LIFT AND DRAG OF 
FLAT SWEPT-BACK WINGS WITH INTERACTING 
LEADING AND TRAILING EDGES 
By Doris Cohen 

SUMMARY 


The problem is considered of a wing with rectilinear plan form swept 
so that both leading and trailing edges lie within their respective Mach 
cones j moreover, the Mach lines from the trailing— edge apex intersect the 
leading edge. Formulas and design charts are presented for the lift in 
such a case, based on approximate formulas for the lift distribution 
developed in NACA TN 1991, 1949. The charts cover a practical range of 
aspect ratios and plan forms of moderate taper, with tips parallel to the 
stream. The leading-edge thrust and drag due to lift are also readily 
calculated from the material presented. Numerical results and an applica- 
tion of the charts are included. 


INTRODUCTION 


Range of Applicability of Formulas 


Problems in linearized supersonic wing theory are characterized 
primarily by the orientation relative to the wing boun dar ies of the Mach 
lines arising from the various points of disturbance defining those 
boundaries. Therefore, even when the plan form of a wing is specified, 
a series of problems arises if the wing is to fly through any consider- 
able range of Mach numbers . 

Paradoxically, the problem of the flow at higher Mach numbers has 
proved more readily amenable to solution than that at lower speeds, 
because at the higher speeds the zones of influence (Mach cones) are 
narrow nnri interference problems are fewer. In the case of the conven- 
tional swept-back wing at an angle of attack, four Mach number ranges may 
be said to have been investigated to a sufficient extent that the lift 
and drag of any specified wing is readily obtainable. The aerodynamic 
characteristics were first computed (reference l) for that range of Mach 
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numhp-rs in which only the center sections of the wing are included in the 
Mach cone from the leading-edge apex; formulas have also been given 
(reference 2) for the next lower range, where the leading-edge Mach line 
cuts the tip; formulas (reference 2 ) and charts (reference 3 ) have been 
given for a third range, where the leading-edge Mach line stands ahead of 
the wing, but the Mach line from the trailing-edge apex is behind the 
wing; and formulas have been given (reference 2) for the fourth case, in 
which the trailing— edge Mach cone includes part of the tip, but not the 
leading edge. 

The last of these cases is illustrated in figure 1(a). In figure 1(b), 
the same wing is shown at a smaller Mach number such that the trailing- 
edge Mach line now intersects the leading edge. In reference *4-, the flow 
in such a case was investigated and formulas for the load distribution 
were derived. In the present report, the load distribution will be inte- 
grated to obtain formulas for the total lift. Formulas for the leading- 
edge thrust and the drag due to lift will also be presented. 

The Mach number range to which the formulas and charts apply is indi- 
cated as a function of taper ratio X, and aspect ratio A for several 
angles of sweepback A in figure 2. The upper-limit curves are determined 
by the condition that the trailing-edge Mach line intersects the wing tip 
at the leading edge. Above these curves, the formulas of reference 2 may 
be used. A lower Mach number range, merging into the transonic and such 
that the flow field of one wing tip extends laterally to include a part of 
the other tip, is defined by the lower— limit curves. The flow pattern in 
such a case is too complicated to be treated by the method of the present 
report. Mach numbers low enough to permit wing— tip interference fall 
within the scope of the approximate theory of reference 5. 

The problem of wing— tip interaction also, of course, sets a lower 
limit on the range of aspect ratios covered. On the other hand, the method 
of reference 4 and the present paper may in theory be applied to the com- 
plete range of aspect ratios to the right of the boundaries shown in 
figure 2. However, the necessary formulas have been developed only up to 
a certain point beyond which the services of high-speed computing machin- 
ery would be required for their evaluation. Both mathematical and physical 
considerations make the application of the method to very highly tapered 
plan forms inadvisable. For cases of moderate taper and a practical range 
of aspect ratios, charts have been computed to facilitate the calculation 
of both lift and drag due to lift. The formulas and charts have been 
developed for wings with tips parallel to the stream only, but the general 
procedure is applicable to wings with raked tips also. 

The symbols used are listed in the appendix. 
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General Outline of the Procedure 


In reference k it was found that, when interaction takes place 
between the flow fields of the leading and trailing edges, the wing plan 
form appears to comprise two principal regions, separated (fig. 1(b)) by 
the Mach line arising at the point of intersection of the trailing— edge 
Mach line and the leading edge. Ahead of this line (region i) the flow 
is most readily described in terms of conical fields. Behind this line 
(region II), the flow is more nearly two-dimensional, and the loading can 
be approximated by the well-known subsonic flat— plate formula, corrected 
in magnitude to give the proper value of the pressure at the leading edge 
of each section. The formulas for correcting the two-dimensional loading 
are given in reference 4. For an infinitely long wing, of course, the 
loading approaches the subsonic flat— plate loading, corrected for the angle 
of yaw according to simple sweep theory. 

The total lift is found using for the loading in region I the conical- 
flow solutions of reference 2 and in region II the formulas given in 
reference From the integrated lift and the leading-edge thrust, the 
drag due to lift can be calculated. The leading-edge thrust in region I is 
the same as that for the triangular wing with the same leading-edge sweep, 
and iias been given by Jones (reference 6), Hayes (reference 7), Robinson 
(cited in reference 8), and others. The disturbance arising at the 
trailing edge affects the leading-edge suction in region II. Formulas are 
derived for calculating the resulting thrust, from which the drag due to 
lift follows. 


FORMULAS FOR LIFT 
Region I 


Conical— flows method .— The term "conical— flows method" is used to 
designate briefly a method of superposition of conical flow fields 1 devel- 
oped by Lagerstrom (reference 9) as an extension of some earlier work by 
Busemann (reference 10) . Thus, the load distribution on a swept— back wing 
is calculated by superposing on the conical loading of an infinite triangu- 
lar plate other conical flow fields of such magnitude and orientation as 
to cancel loading on the portions of the triangular plate lying outside the 
boundaries of the specified swept— back wing. This process, as applied at 
the trailing edge of a swept— back wing, is more fully described in 
reference 2. 


1 A conical field is one in which the velocities are constant along rays 
emanating from a point (the apex of the cone). 
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In figure 1(b) , region I has been further divided into regions la 
and lb by the Mach line from the trailing— edge apex. Region lb is the 
region in which the loading is modified by the addition of the fields used 
to cancel the load behind the trailing edge. In region la, the loading is 
simply that on a triangular wing having the same apex angle as the swept— 
back wing. 

The conical flow fields were originally derived most conveniently as 
velocity, rather than pressure, fields. The loading is obtained from the 
velocity by the linearized relation 

Ap _ Au 
9 " V 


where u is the streamwise component of the perturbation velocity on the 
upper surface, V is the stream velocity, and Ap/q is the coefficient 
of lifting pressure. 

Load distribution in region I.— In terms of the conical variable 
a = Py/x, the loading over the triangular wing is given by 


where 


U A = 




u o = 


mVa 


pE( J 1-hu 2 ) 


( 1 ) 


( 2 ) 


is the streamwise component of the perturbation velocity along the wing 
center line (a=0) and u^ is the same velocity component elsewhere on 
the wing. 

Cancellation of the load in the wake, as described in reference 2, 
is accomplished by superposing first a symmetrical field having a constant- 
load region coincident with the wake to cancel the perturbation velocity 
Uq throughout the wake, and then an infinite number of infinitesimally 
loaded oblique elements to cancel the remaining lift. Figure 3 shows the 
constant— load region (shaded) of a single "oblique" flow field. 

Cancellation of Uo induces on the wing a pressure distribution 
(equation ( 51 ) of reference 2, corrected) proportional to 


(Au)_ = - 




F ( cp) 


K( V 1— mt 2 ) 


(3) 
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where K( *■/ 1— mt 2 ) is the complete elliptic integral of the first kind 

/ 2 

of modulus V 1— mt , and F is the incomplete integral with argument 


cp = sm 


r 1 J 

V 1— mt" 


(M 


Each of the oblique elements induces on the wing the perturbation 
velocity (equation (50) of reference 2) 


(A! 


s -u a _! (l-a)( t a -mt)-(mt-a)(l-t a ) 
,u)„ = a COS X - — 


( 1— nit )(t a -a) 


(5) 


where Uq is the differential of u^ obtained by differentiating with 
respect to a, or 


u„ = 


mu 0 a 


( m 2-a2)3/2 


da 


( 6 ) 


Uncorrected lift in region I .— In order to obtain the total lift in 
region I, 2 we first integrate the uncorrected loading of the triangular 
wing over the entire area, obtaining 


km 


qa p Va 


Co 


f 1+m 
l(l-m)' 


L/ 


cos 


m ' +&2 


m 2 


l-m* 


m(l+a 2 ) 


m — a 2 
l+a 2 






2 2 


(- 




-a2 


HL|- — a 2 


X 


% 


Vm^ 2 — m 2 


3 r 


cos 


_ x 


mi 


(m^.— a^ 


cos 


— l 


m 


\ J 


(7) 


where a 2 is the value of a corresponding to the ray through the point 
of intersection x 2 ,y 2 of the trailing edge and the Mach line separating 
regions I and II. (See figs. 1 and 3-) In terms of the wing plan— form 
parameters , 


a 2 


2m m^. 


1+m+m^-inm^ 


(8) 


formulas which follow are for the complete wingj that is, the left-hand 
side is included. 
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When =m (untapered wing), the second part of equation ( 7 ) becomes 
indeterminate. In this case. 


L n u n „ f 1+m f 1 m 2 +a 2 

v — — ai C ) — - - 0OS 

qa P Va 0 L(l— m) 3 L 1— m 2 m(l+a 2 ) 


y 2 2 

in — a 2 
l+a 2 


+ 


1 

3 



m 

m-a 2 


m+a 2 

m-a 2 



(9) 


Wake corrections .— The wake corrections are to he integrated over 
region lb. Integration of the symmetrical wake correction (equation (3)) 
yields 


where 


Wo = -l6m 2 xx 2 Uq f 
qa PCi+mt) Va \ 


_2 

l-^t 


E ( ") 

K( aA^ 2 ) „ J 


*1 


o 

1-m 


( 10 ) 


( 11 ) 


is the x coordinate of the intersection of the trailing-edge Mach line 
with the leading edge. (See fig. 3*) 

For each oblique element, the reduction in lift is given by 


1 dAL , 

da 

qa da 


2uo 1+m 


3Va 1+a 


(mt-a)(x2-x a ) £ 


( l+mt)(l-a) 


- 1 


( 12 ) 


in which 


x 2 = 


1+m 

1-m 


“t 


i +m t 


(13) 


is the x coordinate of the intersection of the Mach line from x 1 ,y 1 
with the trailing edge (see fig. 3) anl 


x a 


_ 


“tf® 

is the x coordinate of the apex of the element a. 


( 14 ) 
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Lift in region I .- The total lift in region I is then given hy 

fL) ,ht + . L T" 2 ^ da 

\qa/j qa qa qa J Q da 


p ✓ v 

The quantity | --- — ) is plotted against rat in figure 4 for 

m 2 c 0 2 V<l c Vi 

several values of the ratio m/m This latter parameter is the ratio 
of the tangents of the semiapex angles of the leading and trailing edges 
and is constant for any one wing through the Mach number range . 


Eegion II 


Tip-effect problems .— In figure 1(b), region II is shown divided 
into two sections by the Mach line from the leading edge of the tip. 

Such a division is actually an oversimplification of the problem. While 
region Ila can be treated, by the method of reference 4, in such a way 
as automatically to satisfy the Kutta condition at the trailing edge, the 
tip effect, which modifies the lift in region lib, has not been determined 
so as to take this condition into account. In the cancellation method 
used to determine the tip effect, a so-called "primary tip correction" is 
first obtained by superposing conical flows to reduce the lift to zero 
along and outboard of the side edge. A 
further succession of steps (reference 2) 
is required to cancel lift introduced in 
the wake behind region lib, and thereafter 
outboard of the tip within the Mach cone 
from x*,y* (see sketch) and so on. The 
corrections obtained by this procedure 
alternate in sign and become successively 
smaller in magnitude, while increasing in 
mathematical complexity. Since, moreover, 
it is known from experiments that the 
assumed flow in the tip regions is at var- 
iance with the physical flow, it would be 
illogical to attempt any precise evalua- 
tion of these corrections. A simple for- 
mula will therefore be given, following 
the derivation of the primary tip correc- 
tion, for obtaining a fair estimate of the 
"secondary correction"; further correc- 
tions will be neglected. 

Loading in region II .— Except for tip 
losses, the lift in region II may be 
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calculated by integrating out to the tip the approximate load distribution 
given in reference 4. This loading is 


which is merely the subsonic two-dimensional flat- plate loading, applied 
on sections taken normal to the stream, and with a correction factor 
inserted to bring the loading at the leading edge into agreement with that 
derived for the swept-back wing by the conical-flows method. The theoreti- 
cal loading at a subsonic leading edge is, of course, infinite in any case, 
with the infinity entering as the reciprocal square root of the distance 
to the leading edge. The coefficient 0 is a measure of the strength of 
the leading— edge singularity in the loading on the swept-back wing. 

It is most convenient to describe the leading-edge singularity in. 
terms of the coefficient in the perturbation velocity u of (m x-Py)“i/ 2 , 
where the quantity (mx— (3y) is |3 times the spanw ise distance to the 
leading edge. This coefficient is then u-/mx-|3y and the strength of 
the singularity in u is the limit of its value as Py approaches mx 
(the leading edge). The coefficient 0 is a nondimens ional form of the 
same quantity: 


and is a function of x. 

In line with the general procedure of the conical— flows method, the 
coefficient of (mx-py)~ 1//2 corresponding to the triangular-wing 

velocity distribution is found first, and is then corrected to take 
account of the effects of canceling the loading in the wake, by the sym- 
metrical canceling element and by the oblique elements. The resulting 
expression for 0 is 




(17) 


cr(x) = 



a ° dAC 


(18) 


Va, */ c Q 


o 


da 


where 


n / mx 

C A = u o J — 


( 19 ) 


NACA TN 2093 


9 


is the uncorrected (i.e., triangular-wing) coefficient; 


(AC ) 0 = 


-4mu 0 


itm^. 1+m 


K(k)E(k,^)-E(k)F(k,i|r) 


( 20 ) 


with k 


/2m t (l-T Q ) 

'/(l^n t )( To +m t ) ' 


= sin 


—1 


/t^+; 


o +m t 

2t 0 


and t q = 


mx 


x— c. 


is the symmetrical wake correction; and 


dAC 

da = 

da 


-4m u a 



Vl— a*/ (m^.— a)x— m^.c c 


r /(•^- +a )( T a ~ a ) 

Iv (l-a)(m t -a) 


K(k a )E(V+ a ) ~ 


E(k a )F(k a ,i|r a ) 



K( k a )E( k a , \j/ 0 ) 


E(k a )F(k a ,t 0 ) 


(21) 


is the correction for a single oblique wake element, with 


= (mt-ajgg-fflfcCpa = /(m t ^)(l^r a ) 

(m t -a)x-m t c 0 a 7 ( T ft -a ) ( l+mj. ) 


k a = 


(l +m t)(l~ T ft) 

(l-m t )( 1 +T a ) 



The upper limit of integration is 


( 1— m)x— c Q 

0 ^ (l-m)x-m t c 0 


( 22 ) 


Equation (l8) for the coefficient of the leading— edge sigularity is 
mathematically complete for points immediately beyond along the 

leading edge, but if the aspect ratio is very high or the Mach number 
very low, the Mach line from may intersect the leading edge 
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(see sketch) and bring about further 
modifications of the flow on the 
outboard portions of the wing. The 
only effect of this modification on 
the formulas presented herein is the 
introduction of additional terms 
into the expression for cr for 
values of x beyond x 3 , and in 
still more extreme cases , beyond 
x 5 , and so on. Evaluation of these 
terms by the conical-flows method 
involves multiple integration and is 
impractical except with the aid of 
high-speed computing machinery. 

Each successive correction (to a) 
is initially zero and enters with 
zero slope at Xq, zero slope and 
curvature at x 5 and so on, so that 
the three— term expression for a 
given by equation (l8) may be used 
with satisfactory accuracy for some 
distance beyond the last value of 
x (x 3 ) for which it is strictly 
valid. In practice, the third term 
in equation (l8) may also be neg- 
lected for values of x only slightly 
greater than x lf 


Charts giving ^ as a function of X - — ■ ■■■ for the values of 

0 ^ m V m c Q 

the taper— ratio parameter covered in figure k are presented in 


figure 5 as an aid to computing. The curves were computed using equa- 
tion ( 18 ) and are therefore exact only up to x=x 3 (shown by a vertical 
mark in each curve). The points x=x are also indicated (by X f s) as a 
more practical limit to which use of the curves may be extended. (These 
points are off the scale for mt = 0.8 and 0.9 in figure 5(a)). When the 
wings are untapered (m/m^. = 1.0), asymptotes 



derived from simple sweep theory, may be drawn. The relative positions 
of the curves in figure 5(a) and their asymptotes suggest that a further 
extension beyond x 5 will probably not introduce any serious error in 
these cases . 

The curves are for the most part regular enough to permit interpola- 
tion within intervals of 0.2 of mt . However, at m-t = 1.0 the lines 
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diminish to a point on the vertical axis; a curve for m^. =0.9 was 
therefore inserted in the charts for moderate values of m/m^. With 
rn/m-t less than 0.% m^. = 0.9 represents, if x >xj_, such extreme 

taper that the successive points of reflection of the Mach lines 
(x 3 , take place within a very small fraction of a chord— length 

and no useful curve can he drawn. No curves are drawn for mt smaller 
than 0.2 because of the tip— interference limitation discussed in the 
introduction. 

Uncorrected lift, region II .— In order to find the total lift (except 
tip losses) in region II , a~ double integration with respect to x and 
y is performed on equation (l 6). A first integration, with respect to 
y, yields for the indefinite integral 




/ 


m, c — ffl,x+3y 
tot 17 


mx-Py 


(23) 


The values of 0y to be sub- 
stituted as limits in equation (23) 
are indicated in the sketch to the 
right. Along the leading edge, the 
right-hand member of equation (23) 
reduces to zero; along the trailing 
edge it becomes 




Then the total lift in region II 
(on both wing halves), except for tip 
losses, is 
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m 


where 



f 2 = J ( 1+m) ( x-Xi ) f 4 = //m(x— Ps/m) 


f 3 = y(l+n4.)(x2-x) f 5 = 


The indicated integrations may be performed numerically or graphi- 
cally, using values of a(x) taken from the charts. 

Primary tip correction .— A method of approximating the tip— induced 
reduction in loading is described in reference 4. The uncorrected 
loading along the tip section is determined by equation (l6), simplified 
by replacing a(x) by the constant 



(26) 


The assumption is then made that the lift to be canceled outboard of the 
tip is a continuation of this loading along lines parallel to the leading 
edge of the wing. The cancellation is accomplished by means of the "tip 
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solutions” of references 2 and 9> 
flow fields having constant 
pressure over a region hounded hy 
the tip and a ray extending 
outward from a point x c ,s on 
the tip. 

In canceling the cylindri- 
cal field assumed in this case, 
the rays would all he parallel 
to the leading edge. (See 
sketch.) The streanrwise com- 
ponent of velocity on the wing 
due to one such element is 


“c 

— cos 

it 


m+t c +2nrt c 
t c — m 


(27) 


where u c is the constant value 
of the streamwise perturbation 
velocity on the element and 
t c /p is the slope of a ray from 
its apex x c ,s. This velocity 
is multiplied hy 2pV to 
convert to pressure and hy 



dS _ mt 2 [ct-(x c - ^-)]2 

dt c dtc 3(mt-t c ) 2 


to obtain the lift on the differential of area shown in the sketch. Then 
integrating with respect to t c from —1 to 0 gives the total lift induced 
on the wing hy the single canceling element at x c . The result may he 
written 

1 dAL a Suc^^t 2 

dx c = 

qa dx c 3Va( 

where | c , as shown in the sketch, is the distance of x c ,s behind the 
leading edge, divided hy the tip chord c t , or 

x c -(ps/m) 


m +m 2 m \ . . . 2 

. U-U 


(28) 


Co+CPs/mtMPs/m) 


( 29 ) 
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Since along the tip, from equation (l6). 


u(x c ,s) 


= cr. 


Va 


! [3s-m t (x c -c 0 )]c 0 

[ ®tc 0 -( mj.— m)x c ] ( mx c H3s ) 


or 


u(g c ) 

Va 


1— £ , 


\/mX 






(30) 


where X is the taper ratio ct/c 0 > the canceling velocities u<; must 
be given by the differential 


u c . 


Cc d 


1-L 


Va a/eA d| | 


/ n^-m 

l<5 V *t lc 


<U, 


(3D 


Then the primary tip correction to the lift (both tips included) 



H 2 °t 2 / / m+m 2 

P(mt-ffl) VV mt+mt 2 


r 










C 


) 


d *c 


(32) 


which may be integrated by parts to give 


AjL \ — l6cr s m t 2 ct 2 • Jo ./’ / m+m 2 

- qay tip 3P(mt-m) 3 \/x V mt+m^. 2 



(33a) 
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where K and E are the complete elliptic integrals of the modulus 



In the case of an untapered wing, equation (33a) reduces to 


AiL\ _ — 3*c 0 2 \/"m 
qa /tip ~ 2 l 3 ( 1+m ) CTs 


(33b) 


Secondary tip correction.— The 
primary tip effect Just derived is 
actually am overcbrrection, since 
the load distribution of which it 
is the integral does not go to zero 
at the trailing edge. The residual 
lift in the wake is the result of 
superposition of an infinite n umb er 
of conical fields and therefore 
cannot be canceled identically in 
any simple way. Numerical results 
indicate, however, as pointed out 
in the previous reports (references 
2 and 4), that the major part of 
the tip effect arises in the can- 
cellation of the infinite pressure 
at the leading edge. In determin- 
ing the secondary tip correction to 
the loading, therefore, it may be 
as s\imed as an approximation that 
the residual velocity field in the 
wake is conical with respect to the 
leading— edge tip, the value of the 
velocity along any ray t m (see 
sketch) being determined by the value 
at the intersection 
Kutta condition is 
in the wake. 



x b->yb of ^e ray with the trailing edge. Thus the 
satisfied, although some pressure differences remain 


An expression for the correction to the total lift resulting from 
this approximate cancellation of the tip— induced velocities at the 
trailing edge has been given in reference 1 as equation ( 71 ), from which 
may be written 
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( s— y* ) 




-Ji 




d(Au/Va) 


dt 


m 




“t 




( 34 ) 


where x*,y* (see sketch p. 15) is the point of intersection of the tip 
Mach line with the trailing edge, t m is £ times the slope of a ray 
from the leading-edge tip, x-j^y-^ intersection of the ray with 

the trailing edge, Au is the velocity to be canceled at 
and Au* is Au(x*,y*). 

Again experience has shown the integral term in equation ( 34 ), which 
involves very lengthy computing, to be considerably smaller than the term 
in Au*. Since only an approximate tip correction is desired, it will 
suffice, therefore, to calculate 


— ) * _ p- ( s -y *) 2 

qa / H Va v J ; 

tip 

From reference 4 , equation ( 45 ), 


/■ 


l+mt 

% 



( 35 ) 


Au* _ ~Os 
Va mki* 


(3 6 ) 


where |* is the streamwise distance of x*,y* from the leading edge 
expressed as a fraction of the tip chord; that is. 


ct 




m 


It is easily determined that 


3 ( s— y* ) = 




and consequently 


I* = 


l +m t 

i°t(l+m) 
m( l+n^ ) 


( 37 ) 


NACA TN 2093 


17 


With these values, and. Au*/Va from equation (36), equation (35) becomes 

(** h ) % 4 5s. 

\qa A , p 1+nu- 


tip 


P 1+m t a/U +mM 




(38) 


Total tip correction .— The total tip correction, to the degree of 
approximation discussed, is the algebraic sum of equations (33) and (38) 


the tip location. Values of 


PVT 

CTsC-t;2 


and 

\ . 

f m 

and 

( AL 


\q a , 

/tip 


have been plotted in 


figure 6, in a form similar to the chart of (^qa^) (fig* *0 • 


Numerical Example 


As a summary of the method, a sample calculation will be outlined. 
The lift— curve slope Cj^ will be calculated for an untapered wing of 

10— foot chord and 40— foot span, swept back 45° and flying at a Mach 
number of 1.08. Then cot A = 1.0, (3 = 0.4, m = 0.4, mt = 0.4 and, 
from figure 4, 


) = 13-30 c Q 2 = 1330 square feet 

I P 

(it should be noted that ^ = cot a) 

With values of a obtained from figure 5(a), graphical integration of 
equation (25) gives 



— ) = 531 square feet 

qa J 

II 

Since the leading-edge sweep is 45°, the tip is at x = 20, and, from 


figure 5(a), a s =0.635^^ ) = 0*5l8. From figure 6, since A = 1, 


— \ = — 1.69 — — — = —219 square feet 

via ,J p 

tip 
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Combining the three components of the lift, we obtain — = 1612 square 
feet and 



DRAG DUE TO LIFT 


The drag due to lift can be found as the integral around the airfoil 
of the incremental pressure puY times the slope w/v of the surface. 

In the case of a thin flat airfoil, this calculation gives merely the 
lift times the angle of attack over most of the surface. However, as is 
discussed in reference 11 in connection with the two-dimensional wing in 
subsonic flow, an infinite pressure acting on the leading edge results in 
a suction force, which tends to reduce the total drag. When the wing is 
swept behind the Mach lines from its apex, a similar suction can exist in 
supersonic flow, since the pressure on a subsonic leading edge is theoret- 
ically infinite. 


Hayes (references 7 and- 12), Robinson (see reference 8), and others 
have derived the formula for the suction force on a subsonic leading 
edge by ass uming the flow near the leading edge to be essentially two- 
dimensional and applying the results of two-dimensional potential theory. 
The simple result obtained in that manner has been verified for the long 
swept— back wing by application of the somewhat different approach of 
reference 13. 


By the two-dimensional approach, the suction force is found to be 
proportional to the square of the strength of the leading-edge singularity 
in the perturbation velocity u. This quantity has already been discussed 
in connection with the lift in region II, where the strength of the 
leading-edge singularity was defined as the coefficient of the inverse 
square root of mx— Py. This coefficient is C A (equation (19)) on the 

forward part of the wing and Ca + (AC) 0 + / da (see p. 8) 


o 

behind the trailing-edge Mach line. From Robinson’s work, the longitudi- 
nal component of the suction force per unit length in the x direction 
may be expressed, for region I, in terms of C A as follows: 


dT 

dx 


pit 

m 


J 1— m 2 C A 2 


and similarly for the remainder of the wing, with C A replaced by the 
corrected coefficient. 
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For a given wing, the total drag due to lift is the result of sub- 
tracting the total thrust from the product of the lift and the angle of 
attack; that is. 


D 




(39) 


or, in coefficient form. 



where S is the total area of the wing. Writing 


£s 


m 


— dx = £2- */l— m 2 T P 
dx m L Jr, 


_ 2 . 

dx + 


(3s 

m 


we obtain 
|3s 

-2-r m 2i dx 

qarj dx 


= it ,/l-m 2 ^ 




Ao 


dAC 

'° T / da 
o 


C A + (AC) n + / ^ da 


■} 


dx y (Ui) 


Xl 2 + r m a 2 dx 


(42) 


so that 


c D =- 

S Lqa LVVa/ 


xi 2 + 


4c ^ A m 


a 2 dx 


m 


( 43 ) 


Xi 


with the final integration to be performed numerically or graphically. 
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APPLICATION AND DISCUSSION 
Lift— Curve Slope 


The lift-curve slope Cp^ has heen calculated for two families 
of untapered wings , with m constant at 0.2 and 0.4, and with varying 
aspect ratio. The results are plotted against the reduced aspect ratio 
PA in figure 7* The circled pqints indicate the actual cases for which 
calculations were made; the faired curves are only approximations, since 
they ignore slight discontinuities in curvature associated with the onset 
of interaction between the wake and leading edge, and higher-order dis- 
continuities associated with successive reflections of the trailing— edge 
Mach lines. 

Two points corresponding to tapered wings with m = 0.4 and 
m^ =0.6 are also included. In one case (pA = 1.6), and in the case 
of the untapered wing with the same span (PA = 1.2), the trailing-edge 
Mach lines did not intersect the leading— edge, and the values of Cp^ 
were obtained by the formulas of reference 1. It should also be men- 
tioned that the remaining values agree within 2 or 3 percent with values 
calculated entirely by the conical— flows method, the former being slightly 
lower than the latter. 

The curves for the untape red w ings may be seen to be approaching at 
the upper end the value 2itm/ Vl-m 2 given by simple sweep theory. At 
the lower end, the curves should approach the origin along the line 

Cp^ = ^ A given by low— aspect— ratio theory (reference l4). As previously 

mentioned, the present calculations cannot properly be extended below 
PA = 1 because of interference between the flow fields originating at the 
tips. However, two such cases have been included for m = 0.2 because, 
with so much sweep, the wing areas affected are small and the interference 
losses should be negligible. Results in these latter cases may be compared, 
because of the small apex angle relative to the Mach angle, with the results 
of the slender-^wing theory of reference 5. This comparison is shown in 
figure 7, although a discrepancy in plan form lessens the significance of 
the agreement. 3 The results of reference 5 have also been plotted for 
m = 0.4, in which case the assumption of extreme slenderness is no longer 
justified and introduces an appreciable error, (it should be mentioned 
that the asymptote for the slender— wing— theor y curv es is below the value 
given by simple sweep theory by the factor Vl— # ) 

3 The theory of reference 5/ while applicable to any swept— back wing lying 
well within the leading— edge Mach cone, has been worked out only for a 
limited family of plan forms, having straight leading and filleted 
trailing edges, and, consequently, a slight taper. 
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Drag Due to Lift 


In figure 8, the drag— rise factor Cp/Cp, 2 is plotted for the same 
families of wings. Comparison is made with a theoretical minimum for 
slender wings in supersonic flight obtained by E. T. Jones in an unpub- 
lished analysis. Using a method similar to Hayes 1 (reference 12) and 
assuming the wing to be narrow compared with the Mach cone, Jones has 
derived a minimum "wave drag" coefficient 



m 


where A x is the aspect ratio defined in the streamwise, instead of the 
spanwise, direction; that is, if l (numerically equal to x^) is the 
ovei^-all length of the wing, 

A x * z2 / s (45) 


The wave drag is to be added to the "vortex drag, " which is the 
induced drag of subsonic flow, calculated from the spanwise loading. 
Using the minimum induced drag calculated from lifting— line theory gives 
as the minimum supersonic drag— rise factor 

± + _&1 (46) 

Cjf TtA 2 tcA x 


It may be seen that the drag rise of the constant-chord swept-back 
wings is fairly close to the minimum, especially at the lower values of 
m for which equation (46) was derived. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Mar. 15, 1950. 
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APPENDIX 


SYSTEM OF NOTATION 


General 


V 

M 

3 

P 

9 

Ap 

a 

L 


free— stream velocity- 
free— stream Mach number 

Vm 2 — 1 

density of air 
dynamic pressure f -rjpV 2 


pressure difference between upper and lower surfaces, or 
local lift 


"D 


angle of attack, 

radians 

lift 


lift coefficient 

(±) 


\<isy 

drag coefficient 

a) 


Wing Dimensions 


o 

c t 

s 

S 

l 

A 

X 


root chord 
tip chord 
semispan 
wing area 

over-all length of the wing in the streamwise direction 
angle of sweep of the leading edge 
taper ratio (c^./c 0 ) 
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A aspect ratio (4s 2 /s) 


Coordinates 


x,y Cartesian, coordinates in the stream direction and across the 
stream in the plane of the wing 

x & ,y a coordinates of apex of oblique trailing— edge element (see 
equation (l4)) 

x b>7b coordinates of point on trailing edge, within the tip Mach cone 

x c ,s coordinates of point on tip; apex of tip element 

X!,yi coordinates of intersection of trailing-edge Mach cone with 
leading edge (see equation (ll)) 

*■ 2 , 72 . coordinates of intersection of Mach line from Xi,yi with 
trailing edge (see equation (13)) 

x*,y* coordinates of intersection of tip Mach line with trailing edge 

Xrp,s coordinates of intersection of tip and trailing edge 

I* streamwise distance of x*,y* back from leading edge, as a 

fraction of the tip chord (equation (37)) 

| c distance of x c behind leading-edge tip, as a fraction of the 

tip chord (equation (29)) 


In the following, all slopes are measured counterclockwise from a 
line extending downstream from the apex of the wing or of the pertinent 
elementary sector: 


m 


“t 


a 


slope of leading edge 
slope of Mach lines 


0 cot A 


slope of trailing edge 
slope of Mach lines 


slope of ray from the origin 
slope of Mach lines 



slope of ray from trailing— edge apex 
slope of Mach lines 


_X_ 

x-Cq 
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slope of ray from leading-edge tip 
slope of Mach lines 



y-s 


m 


slope of ray from apex of element a 
slope of Mach lines 



slope of ray from apex of tip element 
slope of Mach lines 


y-s 


x— x c 


the value of a corresponding to a trailing-edge element of 
which the apex lies on the Mach forecone of the point at which 
the load is being calculated (equation (22)) 

aU 2 ,y e ) (equation (8)) 

the smallest value of t Q along which cancellation of pressure 
ahead of the leading edge can affect the point at which pressure 

is being calculated. For a point on the leading edge, 

t 0 = -t 0 = “x/(s-c 0 ) 

the smallest values of t a along which cancellation of pressure 
ahead of the leading edge can affect the point at which pressure 

is being calculated. For a point on the leading edge. 


Streamwise Components of Perturbation Velocity 

basic (uncorrected) perturbation velocity as given by solution 
for triangular wing (equation (l)) 

value of u^ at a=0 (equation (2)) 

constant perturbation velocity on canceling (oblique) sector in 
wake (equation (6)) 

constant perturbation velocity on canceling sector outboard 
of tip 

symmetrical trailing— edge correction to (equation (3)) 

dAu 

da, correction to ua due to single oblique trailing-edge 
element (equation (5)) 

value of tip correction to u at the point x*,y* (equation (36)) 


T 


a 


= t a = (mx-)3y a )/(x-x a ) 
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Arbitrary Mathematical Symbols 


C^ value of coefficient of (mx— By) in at the leading 

edge (equation ( 19 )) 

(AC) q decrement in Ca due to reflection of (Au) Q at leading edge 

(equation (20)) 

da decrement in C A due to reflection of (Au)„ at leading edge 

da 

(equation (21)) 

a leading-edge correction coefficient defined by equation ( 17 ) 

a s value of a at leading-edge tip 


Elliptic Integrals 


k modulus of elliptic integral, defined where used (also with 

subscripts) 

cp or >|r argument of elliptic integrals, defined where used (also with 
subscripts ) 


F(k,<p) incomplete elliptic integral of the first kind of modulus k 
and argument 9 


K,K(k) 

E(k,cp) 

E,E(k) 


complete elliptic integral of the first kind of modulus k; 
that is, K = F(k, £) 

incomplete elliptic integral of the second kind 

complete elliptic integral of the second kind of modulus k; 
E(k,£) 
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Figure /. - Plan view of swept -bock wing ot two supersonic 
Mach numbers. 



Figure 2. - Range of applicability of formulas and charts. 
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Figure J- Coordinate System and other Symbols. 
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Figure 4- Chart for the computation of the tiff in region I 
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Figure 5.- Charts for determining <r, the strength of the leading-edge singularity. 
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Figure 5- Continued. 
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Figure 5 - Continued. 
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Figure 5.- Concluded. 
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Figure 6. - Chart for the correction of the tiff for tip effect. 



uu 

VJI 


NACA TN 2093 


NACA-Langley - 5-25-50 - 1075 


7 


UO 

Q\ 



0 I L 

0 / 


A 


2 

/9A 


3 


-i 

* 



Cd 

~C? 


Figure 8. - Drag -rise factor 


for untapered wings 
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